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ABSTRACT 

We present high-resolution L'-band imaging of the inner scattered light struc- 
ture of Class protostar L1527 IRS (IRAS 04368+2557) taken with the Gem- 
ini North telescope. The central point-source like feature seen in Spitzer Space 
Telescope IRAC images is resolved in the Gemini image into a compact bipolar 
structure with a narrow dark lane in the center. Two scattered light lobes are 
extended ~1.8" (200 AU) perpendicular to the direction of the outflow and ~2.5" 
(350 AU) along the outflow axis; the narrow dark lane between the scattered light 
lobes is ~ 0.45" (60 AU) thick. The observations are consistent with our initial 
modeling of a bright inner cavity separated by a dark lane due to extinction along 
the line of sight of the central protostar by the disk (Tobin et al. 2008). The 
bright, compact scattered light might be due to complex inner structure gener- 
ated by the outflow, as suggested in our first paper, or it may more likely be the 
upper layers of the disk forming from infalling matter. 

Subject headings: ISM: individual (L1527) — ISM: jets and outflows — planetary 
systems: protoplanetary disks — stars: formation 



Introduction 



In the earliest stages of star formation, the Class 0/1 phases (lAndre et al.lll993l ; lLada 



19871 ). the newborn protostar is embedded within a dense envelope of gas and dust. In the 
Class phase, little or no emission is generally detected shortward of ~10//m due to high 
extinction toward the central object; however, the bipolar outflows from the central protostar 
and disk carve cavities in the envelopes creating scattered light nebulae visible in the near 
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to mid-infrared (1-8/mi). The morphology of the outflow c avities then resu lt from a ram 
pressure balance between the infalling envelope and outflow (jShu et al.lll99ll ). 



Early mo dels of infalling envelopes assumed spherical collapse (IShul Il977l ) along with 
slow rotation (Ulrich 1976 ; Cassen fc Moosman 1981 ; Terebev et al. 1983). w hich were later 
modified for infall from initially flattened geometries ( jHartmann et al. 1996 ). The rotating 
collapse model provides a natural route for the formation of circumstellar disks within the 
protostellar envelope. However, disks in Class objects have be en difficult to study primarily 
because of the surrounding envelope (e.g. IJorgensen et al.ll2009l ) and insufficient resolution at 
millimeter wavelengths to isolate the protostellar disk. In addition, the small-scale structure 
of outflow cavities in Class protostars has also been elusive due to the envelope extincting 
near-IR scattered l ight e mission, shifting the scattered light emission peak into the mid-IR 
jTobin et al.lhoOTL l2008h . 



One of the nearest low-mass protostars in Taurus (d = 140pc), L1527 IRS (here- 



I system given its edge-on orientat ion ( 



Ohashi et al 



Ohashi et al. 


1997 


Chandler & Richer 


2000 


Robitai 


le et al.ll2007 


Tobin et al.ll2008 


Grama io et al. 


2010T). It has been classified as Class protostar ( 


Chen et al. 


1995 


) but could be a Class 



19971 ). which enhances the amount of 



extinction along the line of sight (jTobin et al.l 120081 ). The moderate resolution IRAC im- 
ages of L1527 from Spitzer reveal bright bipolar cavities extending ~10000 AU in radius 
( IHartmann et al.l 120051 ; iTobin et al.l 120081 ) . The outflow cavities are separated by what ap- 
pears to be a large-scale (~1000 AU) dark lane; however, we ob served a central po int-source 
between the cavities within the dark lane from 3.6 to 8/zm. In ITobin et al.l (120081 hereafter 
Paper I), we were able to explain this by constructing a model assuming that the point-source 
was a compact inner outflow cavity, unresolved by Spitzer, connecting with the larger outer 
cavity. We then proposed high-resolution mid-IR imaging to test this model. 

In this paper, we present high resolution L'-band (3.8/mi) observations of L1527 from 
the Gemini North telescope. This is the highest resolution mid-IR image of the inner enve- 
lope/outflow cavity around a Class protostar. The observations confirm our prediction of 
a compact, bipolar scattered light structure; quite similar to the model constructed in Paper 
I but strongly resembling a disk in scattered light. 



2. Observations and Data Reduction 



We observed L1527 with the 8.1 meter Gemini North telescope at Mauna Kea on 2009 
October 24 and 2009 December 11 using the Near-Infrared Imager (NIRI) in L'-band (3.8/im). 
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The camera was used in f/14 mode which provides a 0.049" pix -1 scale; we only used the 
central 512x512 area of the detector for faster read-out which enabled integrations short 
enough to avoid saturation from sky background. We observed with a standard 5-point 
"dice-5" dither pattern with 5" steps and a nearby star was used for tip-tilt correction and 
guiding. At each position 150 coadded images of 0.2 seconds were taken; a total of 35 images 
were taken yielding 16 minutes of integration during the first observation and 176 images 
totaling 88 minutes of integration during the second observation. However, the first 38 images 
in the second observation were unusable due to variable sky background. Observations of 
the standard star HD 22686 show that the seeing was ~0.3" on both nights. 

The raw data were reduced using the Gemini IRAF0 package. Flat field images were 
generated from sky flats constructed by median combining the dithered images. Each on- 
source frame was sky subtracted using a median sky generated from the on-source frame 
and the two images taken immediately before and after. This method of sky subtraction 
compensates for the rapid variations in sky brightness at L'-band. The sky subtracted images 
were then combined using the imcombine task of IRAF using the offsets given in the header; 
LI 52 7 was quite faint in the individual frames and could not be used for further refinement 
of the offsets. The datasets from both epochs of observation were combined yielding a final 
image with 85 minutes of integration. 



We also include data taken with the Infrared Array Camera (IRAC) (IFazio et al.l 120041 ) 



on the Spitzer Space Telescope. These dat a were taken in 24 Febr uary 2004 as par t of th e 



GTO Taurus survey and were presented in iTobin et al.l (120081 ) and lHartmann et al.l (120051 ) . 



Results 



The images of LI 52 7 at L'-band and 3.6/im are shown in Figure 1. The IRAC image 
clearly shows the point-like structure appearing in the center of the envelope between the 
large scattered light cavities. The Gemini L' image resolves the point-like structure into a 
compact, bipolar scattered light nebula at the center of L1527. The two lobes are separated 
by a narrow dark lane ~0.45" (60 AU) wide, consistent with a circumstellar disk shadow. 
The L' image reveals struc ture similar to that of the Hubble NICMOS images of Class I 
protostars in Taurus from iPadgett et al.l (11999I ). However, this is the first time such an 
image has been captured of a Class protostar. 



X IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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The total extent of the bipolar structure in the direction of the outflow (east-west) is 
'2.5" (350 AU) and -1.8" (250 AU) in width. The extent of the scattered light along the 



outflo w and dark lane is much thicker than a prototypical disk (i.e. HH30, iBurrows et al. 



19961 ) indicating that the disk in L1527 is more vertically extended or that the scattering is 



coming from the inner envelope. The eastern lobe is about twice as bright as the western 
lobe; consistent with eastern side being slightly inclined toward us (—85°) (Paper I). The 
orientation of the disk dark lane indicates that the rotation axis of the system has a position 
angle ~85° east of north. We also see extended, low surface brightness outflow cavity emission 
extending away from the inner envelope in the Gemini image. The cavity is quite narrow 
until ~6" (840 AU) from the center; this narrowness was not well resolved in the Spitzer 
image. Past 6", the cavity expands rapidly to very wide angles. The total flux within a 7"14 
(1000 AU) aperture radius is 6.5 ± 0.6 mJy, consistent with the IRAC flux from Paper I. 



3.1. Models 



We use the Monte Carlo radiative transfer code of IWhitney et al.l (120031 ) to interpret 
the observations. The model components are described in Paper I, but we discuss some 
important aspects here fo r completeness. The envelope density structu re is the standard 
rotating collapse solution (ICassen fc Moosmanlll98ll ; iTerebey et al.lll984l ) in which envelope 
density scales with the infall rate M env with p oc R -3 / 2 outside the centrifugal radius (Rc) 
and p oc R _1//2 inside Rc. Conical outflow cavities are carved out of the envelope, their shape 
is defined by a polynomial, the degree of the polynomial is the shape parameter. We use the 
dust model from Paper I for the envelope which has grains up to 1/im in radius. 

The disk is defined by its initial scale height H at the stellar radius R* and flaring 
power law H oc R 13 , radial density profile p oc R~ a , mass, and outer radius. The vertical 
structure of the disk is assumed to be Gaussian with an initial scale height defined to be in 
hydrostatic equilibrium at the dust destruction radius (T d =1600K). Values of j3 = 9/8, 5/4 
are commonly used, corresponding to isothermal scale heights with T^k oc R~ 3 / 4 and R -1 / 2 . 
The gas and d ust of the dis k are a ssumed to be well mixed. The dust opacities for the disk 
are taken from IWood et al.l (120021 ) which were used to model the SED of HH30. The dust 
is distributed in a quasi-settled manner with grains up to 1mm in the disk midplane while 
smaller grains remain extended. The transition between these dust populations is defined to 
be n# 2 =10 10 cm~ 3 or the first scale height. 

In Paper I, we inferred that the apparent point source seen at 3.6pm was in fact a bipolar 
structure seen in reflected light, with the true central source actually hidden by extinction. To 
create this structure, we constructed a dual-cavity model having a narrow inner cavity with 
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an outer cavity offset by 100 AU (see Figure 10 of Paper I for an illustration). Comparing 
the top and middle rows of Figure 2 shows that our prediction was qualitatively correct. 
We next take advantage of the high spatial resolution of the Gemini data to improve our 
understanding of the inner structure. We required that models reproduce b oth the L' and 



IRAC images and the broadband SED and IRS spectrum (IFurlan et al.ll2008l and Pape r I) 



as it is important to use as many constraints as possible (Paper I, iGramajo et al.ll2010l ). 



We first adjusted the parameters of the dual-cavity model to better reproduce the Gem- 
ini observations by varying the inner and outer cavity shapes and opening angles. As shown 
in Figures 2 and 3, our best-fit dual-cavity model reproduces the observations reasonably 
well (note that the images for the refined dual-cavity model are not shown as they nearly 
identical to the disk model as discussed below). The envelope infall rate for this model was 
10 _5 M Q yr -1 ; the outer cavity was offset from the inner cavity by 85 AU. The disk parame- 
ters are H(10AU)= 1.87 AU with = 1.25, H = 0.0332 R„ M disk = O.O5M and R disk = 25 
AU; a model with /3 = 1.125 was also able to be fit, but with double the initial scale height; 
the full parameter set is listed in Table 1. 

We did not consider a disk model in Paper I because the large-scale dark lane seemed 
to be too thick to be reproduced with a disk. However, using the detailed Gemini image, 
we can describe disk parameters which yield scattering surfaces comparable to that of the 
dual-cavity model. The highly-flared disk in this model had H(10AU)= 1.95 AU (H(190AU) 
= 82 AU), with = 1.27, a = 3.0, H = 0.03 R„ M disk = O.OO5M and R disk = 190 AU. The 
envelope infall rate was lowered to 0.8 x 1O" 5 M yr" 1 in order to allow more scattered light 
from the disk to escape through the envelope. We also used the envelope dust model for 
the disk upper layers rather than the default dust model; this increased the U intensity and 
reduced emission from 10 to 60 /iin (see Table 1 for other parameters) bringing the model 
into closer agreement with the observations. As shown in Figures 2 and 3, this model also 
provides reasonable agreement with the observations. 



4. Discussion 



The dual-cavity model in Paper I was justified based on simulations with a wid e-angle 
outflow and a highly flattened envelope density distribution (IDelamarter et al.l 120001 ) . How- 
ever, it is unclear if such a flattene d envelope is present in L 1527, as there is no evidence for 
this in submillimeter observations ( {Chandler &: Richerll2000l ). In addition, the outflow cavity 
on scales larger that 840 AU from the central object does not easily fit within this framework 
due to the rapid widening of the scattered light nebula beyond 840 AU (§3; Figure 1) This 
is difficult to explain with a simple conical cavity because there is no direct line-of-sight 



-6 - 



from the protostar/disk to these parts of the outflow cavity. Therefore we suggest that the 
morphology of the scattered light nebulosity on large scales is not due directly to th e outflow 
but to the morphology of the envelope and/or ambient medium (ITobin et al.ll2010l ). In any 
case it appears difficult to explain the complex shape of the scattered light structures with 
an simple outflow-envelope interaction. 

Given these difficulties, the disk scenario seems attractive. Images of disks in silhouette 
have previously been obs e rved around T-Tauri stars (Class II sources) and Class I protostars 
(e.g. iPadgett et al.lll999l ; iBurrows et al.lll996l ). The small-scale scattered light morphology 



of L1527 bears a resemblance to that of IRAS 04302+2247 and HH30. The problem for this 
model is that the disk must be highly flared, with the aspect ratio of the disk extinction lane 
relative to the extent of the scattered light surfaces being about a factor of two larger than 
other objects. To reproduce the compact scattered light and large-scale dark lane with a disk, 
it was necessary to adopt a high degree of flaring in the disk (/3=1.27), a steep density profile 
(a=3.0), a large disk radius (Rd=190AU), and a small disk mass (Mrf=O.OO5M ). (Note that 
the mass is simply a formal parameter of the model for scaling the density. Extrapolating 
disk masses from fitting scattered light in the disk upper layers is problematic since they 
contain a small fraction of the total mass.) These obser vations and modeling emphasize 
not o nly that disks form early during protostellar collapse ( IJOrgensen et al.ll2009t IVorobyov 
2009 ) but also tha t disks with large radii are able to form during the Class phase (see also 
Enoch et al.l 120091 ) with radii comparable to Class II disks in Taurus ([Andrews fc Williams 
2007h. 



In general, dust growth and/or depletion in upper disk layers relative to ISM conditions 
is needed to repr oduce the scattered l i ght images of edge-on disk systems and the SEDs 
of T Tauri stars Jp'Alessio et al.l boodJDuchene et all boiol : iFurlan et aD l20Q8h . Thus, it 
may be that the postulated L1527 disk has had much less dust evolution and thus has more 
small dust in upper la yers, making it easier to explain the wide extinction lane observed (e.g. 
D'Alessio et~al1l2006h . 



Loinard et al.l ( 120021 ) directly imaged the disk of L1527 and a binary companion sepa- 



rated from the primary by 0.17" (24 AU) with VLA A=7mm observations. These observations 
are overplotted on the L' image in Figure 4. The observed disk is quite compact - only ~0.3" 
(40 AU) in diameter. Though this system is a close binary, we do not expect the companion 
to affect the modeled parameters significantly, because our re s ults m ostly constrain the outer 
disk. The compactness of the disk measured lLoinard et al.l ( 120021 ) is not inconsistent with 
our disk model since the A=7mm observations will only be sensitive to the densest part of the 
disk and where the grains are largest. Their measured disk mass of O.1M is substantially 
larger than our disk model (which is uncertain as discussed previously). Note that the disk 
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mass is strongly dependent on the assumed dust m ass opacities; lLoinard et al.l ( 120021 ) used 
a dust opacity model based on IPollack et al.l (|1994j) which had K7 mm =4.9xl0~ 4 cm 2 g~ 1 (P. 
D'Alessio private communication)). lAndrews &: Williams! ( 120051 ) used K 8 5 OAtm =0.035 cm 2 g _1 
with k oc A^ 1 which yields K 7mm =0. 00425 cm 2 g~ 1 when extrapo lated. The difference between 



the K7 m m values can be attributed to the IPollack et al. 



(| 19941) dust model only con sidering 



grains up to 1/zm in size while the dust m odel used by lAndrews fc Williams! (12005! ) consid- 
ered grains up to 1mm in si ze. Scaling the lLoinard et al.l (12002! ) A=7mm mass to n-j mm from 



Andrews &: Williams! ( 120051 ) gives Mdi S k=0. 0095 M Q , in agreement with the mean disk mass 
of Taurus. 

From the observations and modeling it is now clear that the small-scale scattered light 
in L1527 is most likely due to a disk. However, the assumed parameters needed to fit the 
disk model (e.g. steep density profile) may only have been necessary to enable our simplistic 
treatment of disk structure to work for L1527. Models which solve for the vertical structure 



of the disk self-consistently (i.e. iD'Alessio et al.lll998! ). given the irradiation from the central 
object have a density structure that is not Gaussian (as assumed in our modeling) but falls 
off less rapidly with z because the temperature increases with vertical height z. Indeed, 
even in the context of an isothermal structure, the Gaussian density distribution becomes a 
poor approximation as the ratio of the vertical height to radial distance becomes z/R — y 1. 
These issues deserve further investigation with physically self-consistent radiative transfer 
modeling of L1527. 



4.1. Conclusions 

We have presented high- resolution L' imaging of the inner envelope of L1527 demon- 
strating that much can be learned of inner envelope structure and newly forming disks with 
high resolution scattered light observations of C lass protos t ars. T he observations strongly 



resemble the dual-cavity model constructed in iTobin et al.l ( 120081 ); however, the observa- 
tions can be equally well reproduced by a large, highly flared disk or a refined dual-cavity 
structure. However, the global scattered light structure is difficult to interpret as purely an 
outflow cavity. We therefore suggest that the bright scattered light structure is most likely 
be a vertically extended disk- in- formation with material falling onto it. 
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Table 1. Model parameters 



Parameter 


Description 


Dual-Cavity Model 


Disk Model 


Parameter Use 


R*(R ) 


Stellar radius 


2.09 


2.09 


fixed (Paper I) 


T»(K) 


Stellar temperature 


4000 


4000 


fixed (Paper I) 


L*(Lq) 


System luminosity 


2.75 


2.75 


fixed (Paper I) 


M,(M Q ) 


Stellar mass 


0.5 


0.5 


fixed (Paper I) 


M disfc (M ) 


Disk mass 


0.05 


0.005 


varied (both models) 


h(10AU) 


Disk scale height at 10AU 


1.87 


1.95 


varied (both models) 


Ho 


Disk scale height at R* 


0.0332 


0.03 


varied (both models) 


a 


Disk radial density exponent 


2.25 


3.0 


varied (both models) 




Disk scale height exponent 


1.25 


1.27 


varied (both models) 


Mdisk(M e yr' 1 ) 


Disk accretion rate 


3.0 xl0~ 7 


3.0 xl0~ 7 


fixed (Paper I) 


^trunc (R.* ) 


Magnetosphere co-rotation radius 


3.0 


3.0 


fixed (Paper I) 


F spot 


Fractional area of accretion hotspot 


0.01 


0.01 


fixed (Paper I) 




Disk inner radius 


14.25 


14.25 


fixed (Paper I) 


^disk,max (AU) 


Disk outer radius 


25 


190 


varied (both models) 


R C (AU) 


Centrifugal radius 


25 


190 


varied (coupled to R d isk,max) 


R-enu ,min ) 


Envelope inner radius 


42.75 


42.75 


fixed (Paper I) 


Renu ,max (AU) 


Envelope outer radius 


15000 


15000 


fixed (Paper I) 


M e „„(M© yr' 1 ) 


Envelope mass infall rate 


1.00 xlO" 5 


0.8 xl0~ 5 


varied slightly (both models) 




Inner cavity shape exponent 


1.3 




varied (dual cavity model) 


bout 


Outer cavity shape exponent 


1.7 


1.5 


varied (both models) 


z otl i(AU) 


Outer cavity offset height 


85 




varied (dual cavity model) 


$ open, in ( ) 


Inner cavity opening angle 


15 




varied (dual cavity model) 


0open,o«t(°) 


Outer cavity opening angle 


15 


20 


varied (both models) 


@inc( ) 


Inclination angle 


85 


85 


fixed (Paper I) 


p c (g cm -3 ) 


Cavity density 








fixed (Paper I) 


Pamb(g cm- 3 ) 


Ambient density 








fixed (Paper I) 
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Fig. 1. — Left: L' image of L1527 showing the full range of the outflow cavity viewed. Notice 
how the outflow cavity rapidly widens about 6" (840 AU) from the protostar. Right: Full 
IRAC 3.6/im image of L1527. The region viewed in the L' observations is outlined. 
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Fig. 2. — Comparison of models to the observations. Top row: Observations of L1527 
from Spitzer and Gemini. Middle Row: Initial model of L1527 from Paper I without any 
knowledge of L1527 at high resolution. Bottom Row: Disk plus single-cavity model. Note 
that the refined dual-cavity model images are not shown as they are nearly identical to the 
disk models. 
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Fig. 3. — Model SEDs for LI 527 with photometry and IRS spectrum from iTobin et al. 



(120081 . and references therein.). Photometry taken with apertures of 71" A (diamonds) and 
T.IA (boxes) (10000 AU and 1000 AU) are plotted. The triangle at 3.8/im is the Gemini 
L' flux within 1000 AU. The model SEDs are plotted for the Disk model (left) and the 
refined dual-cavity model (right) with multiple model apertures of 10000 AU (solid line), 
6000 AU (long-dashed line), and 1000 AU (short-dashed line). The model is clearly deficient 
in flux at long wavelengths; however, external heating is not taken into account and the dust 
temperatures fall below 10K in the outer envelope emitting less in the far-IR and submm. 
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Fig. 4. — L' im age of L1527 with 7m m dust continuum observations (contours) from VLA 
observations by lLoinard et all (120021 ) . Notice that the position angle of the long axis of the 
disk is in the same direction as the disk extinction lane, ~-5° East of North. The second 
7mm point source just east of the disk is the binary companion. 



